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ABSTRACT 



An approximate theoretical analysis is made of a phase 
discriminator to be used as a synchronous detector of a very 
narrow-band PM wave in noise. Experimental work is. reported 
which indicates the superiority of this method of detection 
over a conventional PM discriminator. 

A comparison is made of the effect upon limiter opera- 
ting point of two different methods of automatic gain control 
A simple diode detector is shown to be superior to a "carrier 
only" method when the carrier-to-noise ratio varies over a 
large range. 
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that provided by Mr. Arthur E. Noyes and Mr. Homer A. 
Lasitter. The continued assistance and encouragement of 
Professors C. P. Klamm , Jr. and P. E. Cooper of the U. S. 
Naval Postgraduate School is sincerely appreciated. The in- 
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investigation. 
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TABLE OF SYMBOLS AML ABBREVIATIONS 



•^subscript 


Amplifier gain at frequency indicated by sub- 
script 


AGO 


Automatic gain control 


cC 


Arbitrary constant 


Subscript 


Abbreviated form for E g sin £ ^ ( ' t su ^ ) )J 


B 


Bandwidth in cycles per second 


Subscript 


Modulation index at point indicated by sub- 
script 


5(f) 


Lirac Lelta Function 


e a 


L.C. self-bias voltage 


e r,s 


Amplitude of sinusoidal wave at point indi- 
cated by subscript 


Subscript 


Instantaneous voltage at point indicated by 
subscript 


f 


Frequency; function of variables in paren- 
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^subscript 


Frequency indicated by subscript. (f Q = 30 Me; 
f^ = 100 Kc in experimental work; f^ = 491 Kc 
in theoretical work) 


&C J 


Fourier transform of term in brackets 


B( ju) 


Fourier transform of a current vs. voltage 
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G(f), G (co) 


Power spectral density 
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n( t) 


Noise voltage 


^superscript 
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script 


<p(r) 
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V'(t) 


Autocorrelation function in I.F. 
Autocorrelation function of output 
Modulating voltage 
Characteristic function 


r t 
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Mean-square voltage in I.F. amplifier 
Normalized autocorrelation function 


R 




R x 


Grid-cathode internal resistance 


S 1 


Carrier-to-noise ratio (in I.F. amplifier) 


1 

s i 


Carrier-to-noise ratio after bandwidth re- 
duction 


S 2 

cr z 


Signal-to-noise ratio (post-detection) 
Mean-square noise voltage 


T(s 2 ) 


Abbreviated expression for a longer function 
of S 


r 


Time interval between t-^ and t^ 


w 

subscript 


Probability distribution indicated by sub- 
script 


x(t), y(t) 


Quadrature components of n( t) . Values at 
specific times indicated by subscripts, e.g., 
x(t 2 ) = x 2 




Dummy variable used for convenience 


vV^VvVvW 


Time average of the quantity beneath wavy 
line 



Ensemble average of the quantity beneath 
solid line 



Time or ensemble average of the term follow- 
ing 

Approaches 

Equal by definition 



1. Introduction. 



Cubic Corporation of San Diego, California, produces a 
Distance Measuring Equipment of the phase comparison type. 
When used for missile tracking the ground station transmits 
a carrier in the 200-400 Me frequency range, frequency modu- 
lated by five lower frequencies. These frequencies are 
selected to permit unambiguous range determination out to 
one-half the wavelength of the lowest frequency. The mis- 
sile being tracked carries a transponder which receives the 
signal from the ground, extracts the modulating frequencies 
and remodulates a carrier of a different frequency which is 
transmitted to the ground. 

The accuracy of the entire system depends upon the 
phase shift experienced by the modulating wave as it passes 
through the transponder. A constant phase shift can be 
calibrated out by simply adjusting the range indicating 
servo at the ground station to a known range. Thus it is 
desirable to use on I.P. strip which exhibits a linear phase 
shift vs. frequency characteristic over a large range of in- 
put signal level. Van Voorhis has shown that synchronously- 
tuned stages are very tolerant to changes in tube capacity 
and their simplicity makes them desirable for this purpose 
(23). Cubic's I.P. strips have six synchronously-tuned 
stages. The Q of plate tank circuits is deliberately made 
low (about 5.25) to decrease the steepness of the phase 
characteristic at midband and improve linearity at higher 

* Numbers in parentheses refer to bibliography. 
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and lower frequencies; stated differently, the overall band- 
width of the I.F. strip is established by adjusting Q. 

To insure operation within the relatively linear por- 
tion of the phase vs. frequency characteristic, Cubic has 
considered it necessary to decrease the effective modulation 
index to a value sufficiently small that only the 

first-order side-bands of the highest modulating frequency 
are significant. The method of negative feedback developed 
by Chaffee (4» 7, 8) is used to reduce to approxi- 

mately.. 02. The feedback voltage is actually a portion of 
the re-transmitted signal which is offset from the incoming 
signal by the I.F. frequency of 30 He. 

The I.F. strip is followed by a limiter and Foster- 
Seeley discriminator. It has been Cubic’s experience that 
the modulating signal power out of the discriminator must 
exceed noise power by about 10 db to permit accurate deter- 
mination of range. Since this deterioration of signal-to- 
noise ratio marks the limit of useful range, consideration 
was given to methods of improving performance. Several 
suggested changes were collected in an unpublished proposal 
submitted to the U. S. Air Force. Fig. (1-1) is a block 
diagram of the proposed transponder. It is the purpose of 
this paper to investigate certain of these suggested 
changes . 
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2. Characteristics of Random Noise. 

In recent years much has been written about the detec- 
tion of signals in noise. A rigorous treatment of any such 
problem requires the use of probability, since typical re- 
ceiver noise is a random process. 

This process cannot be written as an explicit function 
of time, for even if its exact values over a long interval 
of time are known, future values cannot be predicted pre- 
cisely. It is possible, however, to determine the probabil- 
ity that the process will assume a certain value at any in- 
stant. This determination may be made from the observed 
values taken on by the process over an infinitely long in- 
terval of time, or by calculations when the characteristics 
of the primary source are known. Following the notation of 
Lawson and Uhlenbeck (16), this probability will be called 
the first probability distribution and is written W.(y,t), 
where y(t) is the random process and t is any instant of 

V 

time. The assumption is made that the process is station- 
ary, that is, any initial transients have decayed to negli- 
gibly small values and the process will display the same 
statistical properties over any (equal) long intervals of 
time . 



In a similar way an expression may be written for the 
probability that y(t) will lie in the range y^q-^y^ at 
time t^ and in the range y^ -f at a later time t . 

This probability is called the second probability distribu- 
tion and is written W (y , t ; y , t ) or more simply 

2 11 2 2 
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W 2 (y 1 ,y 2 ; X ), where X - t g 




probability distributions may be written, each describing 
the process in more detail. As Lawson and Uhlenbeck have 
pointed out on p. 35 of their book, this series of distri- 
butions completely defines the process. 

A convenient concept in this work is that of ensemble 
and time averages (16). If the random process y(t) is re- 
corded for a very long interval of time and then this re- 
cording cut into a large number (N) of intervals of length 
Q, this ensemble of recordings may be considered the simul- 
taneous outputs of N identical receivers. Now the average 
value of a function of y at a time t may be written 



where the bar indicates an ensemble average. This average 
is, in general, a function of jt. If the average over time 



wavy line indicates a time average. If the averaging is 
carried out in the reverse prder the result is expressed as 



finite time average of y(t) is obtained in either case. If 




(2-1) 



is performed the result is expressed as -£■ , where the 



-f- . The two methods yield the same answer since the in 



the process is stationary, since 'f(tf) is 

independent of time and '£('%■) is same for every member 
of the ensemble. 
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Of special interest is the average value of y(t)y( t-f-lr) 
Denoting this average by ^(?) , it may be written 

= jr*) ?(*''*■) = • 

is called the autocorrelation function of y(t), 
and for a stationary process becomes 

^ r) = I j ^ h r ) 



-00 



= M 






(2-2) 



% 

0 V W„. 

It is seen that r. mean square value of y(t). 

If the average value of y is zero, Cp(0)-= O' =■ variance. 

The importance of the autocorrelation function lies in 
the fact that $(''?} and G(<^) , the power spectral density, 
are the Fourier transforms of each other (24): 

C/O oo 

Sto? r 






'-(70 

oo 



G(/)ct«.iTr-fr df 



0 



~gR> = \Jn) e JUJr Jr 



(2-3) 



t-J—oo 



G(-f) - y 



(p('l') c*<l ZTf-fr d*? 0 



o -oo 
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In a large class of problems, the probability density 
functions for random noise voltages are relatively simple 



mathematical expressions. Lawson and Uhlenbeck (16) and 
Rice (19) are among those who have shown that typical re- 
ceiver noise sources, such as shot effect, thermal noise, 
induced grid noise, etc., are characterized by a normal or 
Gaussian probability density function of the form 



Chapter 4 of Lawson and Uhlenbeck' s book and Part II 
of Rice's paper contain extensive references to the litera- 
ture. Knudtzon (15) has published experimental verification 
for a variety of primary sources by measuring the autocorre- 
lation function using sampling techniques and then computing 
the power spectrum with an electronic differential analyzer. 

Two different representations of the random noise vol- 
tage will be found useful here. In the first form the vol- 
tage is called n(t), with a first probability distribution 



In the second form the noise voltage is developed in a 
Fourier series occupying a narrow band around the I.F. fre- 



ft (£) - X (£) C**-' Z ?f-f 0 iz ~b ZTT-fc^ ) ( 2 " 6 ) 





quency -f~ 0 , such that (16) 
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where 



Oo 




^=-00 




^--00 



and -/^ is measured from 4- 0 . The a^. and b^ are independent 
random variables and x and y are orthogonal functions such 
that 

X l (t) = tfU) = o' 2 " j 

Y{£) -t f i) = O -far 4// 'l' > ( 2 - 7 ) 



Hence , 

K ( *, V) = 

and 



,2. n <T Z 



£ 



2 2. 

xllzJZ. 

Z cr*- 



e Y[ z^o-f) 

~2f\x,* z t 



where 



N <^r) 

P (r)= 77 ^ 



cr 2 



( 2 - 8 ) 



(2-9) 
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3. Effect of a Linear Device Upon Random Noise. 

It is seen from Eq. (2-3) above that if G{oj) is a 
constant for all values of , then 



and the correlation between successive values of y is zero. 

* 

This limiting case is called the "purely random process" 

(11) and the associated spectrum a "white spectrum". When 
such noise is modified by a linear device so that the output 
spectrum has the shape of the power frequency response of 
the device, the output voltage is still Gaussian (16) but 
the autocorrelation function is no longer discontinuous. 





< ' /?«=' r € 
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4. Effect of Non-Linear Devices Upon Random Noise. 

Lawson and Uhlenbeck (16), Rice (19) > Middleton (17> 

18) and many others have published results of investigations 

; 

of several types of non-linear devices under a -variety of 
signal and noise conditions. 

There are two generally accepted methods available for 
the study of the output of non-linear devices. In both of 
them the output power spectrum, at least in the range of 
interest, is determined. The first method, called by Rice 
the "direct method" and attributed by him to W. R. Bennett, 
F. C. Williams, J. R. Ragazzini and others, is based upon 
the known theory of modulation in non-linear devices. In 
this method a typical modulation product is computed and 
the output power spectrum then obtained by summing all such 
products . 

The second method, attributed (19) to D. 0. North, 

J. H. Van Vleck and D. Middleton, consists of first finding 
the autocorrelation function of the output, then taking the 
Fourier transform to give the power spectrum of the output. 
In this method use is frequently made of the characteristic 
function to more simply evaluate the autocorrelation func- 
tion. The characteristic function is widely used by Lawson 
and Uhlenbeck throughout their book. Rice also uses the 
characteristic function approach in many cases because of 
its ready application to many devices expressible by the in- 
tegral form 
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( 4 - 1 ) 



where F( ju) is the Fourier transform of the I vs. V. charac- 
teristic of the device. 

In many problems one of the methods may prove simpler 
mathematically than the other. 
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5. Analysis of a Phase Discriminator. 



One of the changes suggested by Cubic is to replace the 
frequency discriminator with a phase discriminator which will 
be provided with a locally generated reference voltage in 
phase quadrature with the I.P. carrier wave. Cubic intends 
to supply this reference voltage by using an adaptation of the 
"active filter" designed by Jensen and McGeogh (14) to filter 
a portion of the I.P. signal which will then be shifted in 
phase by 90°. In this paper it will be assumed that the refer- 
ence voltage is available. 

A sketch of the proposed phase discriminator is shown 
below in Pig. (5-1). 




Pig. (5-1) 



/ 
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The voltages applied between the plates of V-l and V-2, 



respectively, and R.F. ground are 



e. - £. 



' f? ~E ~h E $ «y u) c jt -h^/tTl -h H (•£) 



^2. — Eg> C^-o.- u) 0 ^ cO c £ +. ijs(£) j — n ( E) 

or, since ^ is small, 

0 f (^E#, £ E 5 Ao^tp-f- Xj(3<^u) 0 3 -f-(E s ~t o) ti jo 
6^ ~ (e r - E s fi-x) ^\E —( £ s u) 0 ± 

where n(t) hab been expressed as in Eq. (2-6). 

Considering the two diodes to be ideal linear rectifiers, 
the output voltage will be a high frequency wave with an en- 
velope given by 



(5-D 



e o (Eg ~h E^aEl, ty + yf-h [E + y-f 



^ j (Ef- F s 4^u(/--x) z -h(t 5 £-yf ' 



(5-2) 



Expanding the expressions under the radicals and sub- 

/ / 



tracting .z. 

e c - le s ^+zx -h 






- (Ej-f"#) 

3 

4- (hlA 

/6 



i r 



/ 



(Eft ~£E$ ^ -h xf (j; -Xy~ 

i __ / 

(E R £-E s ^ t l'+*f 



+ * 
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Terms after the third are neglected since the condition 



will be imposed E = 5E_ . It is recognized that values of 
JX and | ^-j equal to or greater than E r will occur, but 
within the range of our interest the probability of occur- 
ence will be very small. Hence 



e a = zf 5 ^ -h zx 

z 



Es^i'+X £f-E s ^4-X 



.(5-3) 



The binomial expansion is used again to expand the ex- 
pression within square brackets. 






— — £5 'kZt. iztZ j- 

% £ 

/ . E s (p-h x 

7“ 



(t-$ x) 



E 3 



~h 









-h 



(EsjZeZIi >7 



£3 



+ 



-f Z- ( £s X-Z-is t/' -f~ xj 



£ ^ 

^x 

where terms containing E r to the fourth and higher powers 
in the denominator are neglected. Thus, 

"J£*~ 

2 



e 0 ±Z£ s ^,.f+-zx - 



■LZEs^^-f-ZY -E s \—£<I'+ £+ 



/? Z 






+ ZX% +. + X 

£s Ex 
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Z.'A-y 






' /?* 



£< R 2 ' 



f 2 k 



x y 

Rl K 



/ - 



A* 



£ s R 2 - 



(5-4) 



where R 



— £1 



A 



Now denoting by subscripts the quantities at t^ and t^» 

e.g. , x( t ) - x , y( t ) ^ y , etc . , and letting OC a E Aorcp , 
1 1 2 2 8 ' 

the autocorrelation function of the output voltage may be 
written 

(f> o ( R) — <3 0 {£.) <3 0 t£ -h'V)- 



r 



cov<A z x , 



c 



. l&k 



1 - j irJW'Af 14 ?- —Ez'- - *,»? ■} 

. n t <J t t H L E S R* E * R Z £ i R zj 

- iL-2xJ 1-jJ- - I*'-.- 1 - A 3q (5-5) 

if! >1 *\ £sR l 4^ 4V 



Twenty terms (of the 36) are zero because 

— ■ . . — - ■ 

= = *< = *z = y, = '?x - o . 

Thus , 

4>o (V = * - 2 -gf^£ - $ + V v *//- Z 



- ^- X, Xz fe r// - -U- £ h-R\+ t ^ **■ ^ 

£- s/? z L RJ R 7L L Rtl E* R^ 



j ^2^ R ./~j_l/+ R x i *2. &/ y-z 

^ e;v ~rrw~~ 



(continued on following page) 
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^ y, ¥*. */ 
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-h 
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-h 
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It is known that 










xp - xp 


- -- 


fr* = 




Also 


f 










x,x, = 




\ 


The 


noise in the I.P. 


amplifier is 


taken to 


have a rectangU' 



lar power spectrum of width B. Hence, 

o=> 





jZTT-fr ' 

df = 



_ oo 




CT 



ZTrjTB 



' jzrr-fr 

£ 



*/z 



- (T 
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]-e/z 



7Tlr 3 



and 
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Since x and y are independent, their joint probability 
distribution function may be written as a product , e.g., 

j 



^ f X / > <f< ) ^ ) &*> ( ) = 



) ZTT(J' Z '/^ 



‘ 2_Tf (T* /T-p* 









Thus the two terms containing X, X^. ^ and X ( X^ ^ may be 
written respectively X ( X^ 0 an< ^ and are 

both zero. By the same reasoning X,Vx - (T* ^(T) 

ma X,X4, ft, = cr^p^r). 

Calculation of the ensemble averages ^ ^ » 



and ^ is most easily accomplished using the character- 

istic function method of Appendix (A) . Thus 

>/ z ) = £ i + Z f ^ + 



and 



- — < j> < /Y^y , ^4) 

12 (7 x ^ t 



dZ, d z Z z 



z^o 

u^o 



- " (^/ p A " 1 1 ^/ + o 

Z t r .Q 

= o 



Similarly, 



z. _ 

ch <fx “ ^ 
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Thus , 






s*AA 



y^p (///-+ zp*) 



e;/? + e* a 



+ «A„r,-if- yrtp_n + W; . rfaz f) (H . 6) 
+ a' d E;el! *j e*r* ITrv ' 



Before calculating the assumption will "be made 

that only one modulating frequency (AJ ) is present and that 
the effective modulation index is f% . Using the approxima- 
tion 

xl/n. ^ = A — p| ±7 

the calculation is made simpler and a direct comparison may 
be made with the experimental results of Section (6). Thus, 
■2- ^Z. 



^ ^ E S 









7 



|S 7(7 y Ac.-rT-' C4 -) f ( A -h 'X J 



= El ~zPtt ^ ^ r ) 

Jq 

1 ( W 

- — — ) Lu\c/x 

271 \ J* JT 



z 






and 






^ rt- 
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The power spectrum of the output is found by taking the 
Fourier transform of this expression. Using the results of 
Appendix (B) , and writing the spectrum in positive frequen- 
cies only, 





O < f < 4-, 

< f < e-f ( 



e/se iv he re 





o< f <i 

4 - < + 

else^/iere , (5-8) 
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6. Comparison of Theoretical and Experimental Results for 
the Phase Discriminator. 

In the work at Cubic an accurate comparison was desired 
of the relative performance of the frequency discriminator in 
present use and the phase discriminator. Accordingly, labora- 
tory tests were conducted using an actual transponder I.P. 
strip and frequency discriminator. The latter was converted 
to a phase discriminator by disconnecting the end of C-j_ marked 
S (see sketch below) from the primary of the transformer and 
feeding the reference voltage through C^. The reference volt- 
age entered through a shielded cable, the lead of C^ being 
the only unshielded path in both circuit configurations. 



S(£) ~h n (£) 







Pig. (6-1) 



The complete laboratory set-up is shown in Pig. (6-2). The 
noise voltage n(t) was tube noise and thermal noise from the 
I.P. stages; carrier- to-noise (21) ratio (S^) was decreased 
by reducing the signal level. 
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Pig. (6-2) 






(S ) was measured using the two discriminators in succession. 



All measurements were repeated on different days. Pig. (6-3) 
shows photographs of the oscilloscope presentations for vari- 
ous levels of S^. 

Eq. (5-8) may "be used to determine S 2 » • In order to 
compare with the laboratory set-up, let 

-P. = too Kc 

e, = .27 

6 - ZMc 

2 y - Kc 

where Z'ti is the effective noise bandwidth of the post de- 
tection filter and is calculated from its measured 3 db 
bandwidth using the table on p. 177 of Threshold Signals . 

This table also indicates that for the six single-tuned 
stages of the I.F. amplifier, the noise bandwidth closely 
approaches the 3 db bandwidth. 

^ r 

Since S, = and -f - 0 O5D Eq. (5-8) may be used 



2 



to give 




<2 A/o/sg X'/v Output" 
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Pig. (6-3) 
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Curves of vs. for R = 10 and R = 5 are computed 
from Eq. (6-1) and are shown in Pig. (6-4). If the addi- 
tional condition is imposed that 

El « £% » o-*-, 

then the asymptotic condition exists wherein 

-5^= K5, (6-2) 

where K depends upon ^ . This curve is also shown in Pig. 
(6-4). 

Fig. (6-5) is a plot of experimental data showing the 
effect upon of varying the ratio R. 

Pig. (6-6) is a comparison between theoretical and ex- 
perimental results when the conditions of Eq. (6-2) are met. 
Pig. (6-7) is a comparison of phase- and frequency- 
i discriminator performance for a smaller modulation index, 

P = .125* The theoretical result predicted by Eq. (6-2) is 
also shown. 
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7. Automatic Gain Control (AGC). 

The present transponder circuitry employs a limiter 
stage and AGC to the first four I.F. stages. The AGC detec- 
tor is an unbiased crystal diode. Both AGC and limiting are 
considered necessary because of the large range of input 
carrier level. The prospect of increasing range even further 
by using the phase discriminator has aroused interest in an 
evaluation of different types of AGC. 

Of primary interest in this evaluation are the relative 
effects upon phase shift through the I.F. amplifier and lim- 
iter as the carrier level varies over the range. The change 
in phase shift through these components is attributed to 
variations in the input capacity with bias voltage. A change 
in bias voltage alters the gain of the stage (which in turn 
changes the input capacity due to the Killer effect) and the 
direct plate current (which changes the grid-cathode imped- 
ance due to its effect on the size of the virtual cathode) . 
The present I.F. amplifier is considered sufficiently insen- 
sitive to these variations that it can be used with the 
present AGC or either of two proposed types of AGC. Accord- 
ingly, the manner in which limiter operating point varies 
with the strength and composition of the signal at its grid 
will be considered. 

The input circuit of the limiter stage is shown in Fig. 
(7-1). If the path from G to ground is regarded as an ideal 
linear rectifier in series with a fixed resistance R , an 

A 

equivalent circuit may be drawn as in Fig. (7-2). 
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5702WA 

z--_: (similar to 6AK5) 



Pig. (7-1) — Limiter input circuit 



HI— 





Pig. (7-2) — Equivalent circuit 

Since the grid path will conduct in only one direction 
and the capacitor will pass no direct current, any direct 
current flowing through the crystal must return to G through 
the 3K resistor. The d.c. grid bias is thus 

- Rc^ Jx d . c< • (7-1) 

The presence of creates a biased linear rectifier. 
Let the input voltage be 

5(£) -h n(£) - +- n(£) ( 7 _ 2 ) 
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where the signal portion is considered unmodulated to make 
calculations more simple. Since the actual signal is fre- 
quency-modulated and a time average will be taken to obtain 
the d.c. component, the result obtained by using Eq. (7-2) 
will be correct. 

The noise voltage is taken to be Gaussian with a first 
probability distribution 




<T*fzn 




(7-3) 



Current flows through the diode when S(t)-f n(t) - E A >0 
where is the d.c. bias voltage. Replacing S ( t)-f- n( t)- E^ 
by the random variable Z and using the relation (9) 

[*+*>- saQ; (7 _ 4: 



where f(n) is the first probability distribution of n, the 
first probability distribution of Z is obtained, 




[2 + £*-S(tj/ Z 

__L_ € 2 .^ 

<T '\fzrn 



(7-5) 



Now -Z 1 — — where the ensemble average 



* Xd 

is taken over the region S( t)-j- n( t) -E A > 0, or 



T I 

- L x . = -vr zr 



•d.c. 



(7-6) 



over all positive Z. The ensemble average is taken first 
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z = 
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s(£)-Ea 



sU. 



o~ izn 1 ^ zcr ° // '^ c * 

V^) 



The first integral on the right is integrated directly to 
give 

lsw-e a ; z 

cr 




^[ZJT 




and using the relation (12) 



CO 




€ ^ = 



V T T 

T" 



the result is 




(7-7) 



Prom Eq. (7—1) and (7?-6) a conditional equation can be 



written for E^. 
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» (7-8) 

cr v * (T / 
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It does not appear that Eq. (7-8) can be expressed in closed 
form. Bennett (1), on p. 170 of his paper, implied that nu- 
merical methods are necessary to solve this problem. However, 
for a given value of carrier- to-noise ratio, the numerical in- 
tegration required to determine the time average is not too 
tedious. Values of 1^ c for -10db<^S^ <y<-10db are of primary 
interest, since Middleton (17) and Rice (19) have written ex- 
pressions for, respectively, the no-signal case and very large 
carrier- to-noise ratios. In fact, with S(t) = 0, Eq. (7-7) re- 
duces to Eq. 4-3 > P* 785 of Middleton’s paper. It will be 
shown later that for large S^, results obtained from Eq. (7-8) 
are very nearly the same as those given by Rice on p. 123 of 
his paper (Vol. 24). 

An example will be used to illustrate the method of nu- 
merical integration. Starting with Eq. (7-8) and choosing 

S. = ; - / , the result is 

I 2cr ,< - 








e z dt} os) 



v o'O 

where use is made of the fact that - Q . 
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In the actual circuit R is 3KX1 . An arbitrary value of 

O 

1K1L is assigned to R . This will affect the numerical re- 
suits but not the general conclusions to be reached. A value 
is now assigned to and the two terms to be averaged 

are plotted on graph paper from ^ to -h in steps of 

0.2 radians. The average value of each term is computed from 
the number of squares under the curve and the right side of 
the equation is thus reduced to a fraction of E^. The process 
is repeated for different values of until the right side 

is just E a . 

Values for E^ were computed in this manner for carrier- 
to-noise ratios of zero, six, and ten decibels. In the tabu- 
lation below, E a is expressed as a fraction of P^, defined by 




( 7 - 10 ) 



and equal to the r.m.s. power which would be dissipated by the 



input voltage if 


it were impressed upon a one-ohm resistor: 


S x (db) 


e a 




Less than -10 


•552 P t 




0 


.565 P t 




6 


.600 P t 




10 


(.601 P t 


(Numerical Integration) 




i- 607 p; 


(Rice's Approximation) 



The conditions for Rice's approximation are not fully met, 
since he required Eg- 0“ and for S]_ = lOdb, Eg - E^ =3<3~. 

An inspection of this tabulation reveals that the input 
voltage may vary from "pure carrier" to "pure noise" with an 
accompanying variation in limiter self-bias of less than 10 
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percent, provided that the mean square power input to the lim- 
iter is kept constant. 

Both of the proposed AGC systems are based upon circuitry 

which would present an extremely narrow bandwidth to noise. 

The system in Fig. (1-1) utilizes a crystal oscillator to 

translate downward in frequency a portion of the I.F. voltage. 

At this lower frequency (100 K c ) the bandwidth may be reduced 

conveniently to 15 K Q . The carrier- to-noise ratio after the 

1 

bandwidth reduction (S^) is then related to S-|_ by 

/S * /o 3 -5, 



or, 



5/ =. / 33.3 -S', 

S/(c/6)~ S,{db) +21.7. c/6 . 



Thus the AGC detector is only slightly affected by noise until 

the carrier is greatly reduced below noise. 

\ 

The other system is also a "carrier-only" type. Here the 
reference voltage would be shifted in phase by 90° 

and multiplied by a portion of the I.F. voltage. The output 
circuit of the multiplier would have a very low cut-off fre- 
quency, hence a low equivalent noise bandwidth and the d.c. 
output would be proportional to the carrier level in the I.F. 
voltage wave. 

The simple diode which acts as the detector in the pres- 
ent AGC system has been investigated quite thoroughly by Rice 
(19) • His Eq. (4-1-14) shows that for a quadratic detector 

I d .C ~ + fJ (7 - n) 

u > 

where c£ = constant. 
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(7-12) 



Using the notation of Eq. (7-10) this may be written 

^d.c. = « Pz- 

Similarly for a linear rectifier Rice obtained 



T = £ 

tfjjf 



JL 

z. 



(7-13) 



where I and I are Bessel functions of imaginary argument. 

£-2- 

Again using Eq. (7-10) and the relation 5-= — s Eq. (7-13) 

1 Z<r^~ 1 



may be written 






vT+s, 



(7-14) 



In Appendix C it is shown that the quantity within square 
brackets is very nearly unity over the entire range of 
from 0 to c> o . Hence, for the linear detector 

Zd.c. == X / w 

If ideal AGO action is assumed, it is apparent that the 
crystal diode detector tends to maintain the total power out- 
put from the I.F. amplifier constant. The "carrier-only" AGO 

tends to keep E g constant. From Eq. (7-10), by substituting 
/r* 2. 

-2- -=. c~ the result is ' 






K-iX' + i). 



(7-15) 



For S-^<1, decreasing requires that increase at a 

’ -2 



rate proportional to S. 



Cubic anticipates that eventually 



the equipment will operate down in the region where S-^ ap- 
proaches -20db. It is apparent that in a system which has 
appreciable gain for large carrier-to-noise ratios, the re- 
quirement of Eq. (7-15) cannot be met with a practical 
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amplifier for 1. Jaffe and Rechtin (13) have pointed 

out this inapplicability of "carrier-only" AGO to systems re- 
quiring a large dynamic range. 

There is another aspect of the AGO problem worthy of con- 
sideration. Prom Eq. (5-8) it is seen that the signal power 
out of the phase discriminator is proportional to |3 ( . The 
phase discriminator is followed by a "compensating amplifier" 
of fixed gain which amplifies the five modulating frequencies 
before they are used to phase modulate a carrier. This car- 
rier is then multiplied up in frequency until it is 30 Me 
below the incoming carrier frequency, a process which also 
multiplies the modulation index by the same factor (22). The 
modulation index out of the phase modulator at any given modu- 
lating frequency is proportional to the voltage at that fre- 
quency (3)* Hence, the modulation index of the retransmitted 
wave depends upon . But 



moauxaxion index ax r-^ ox rexransmixxea signal. 
Thus, if the circuit employs an AGO system which tends 
to maintain constant, a reduction in S-j_ will result in a 



^ ,OZ • If the phase modulator is assumed to exhibit a 
linear relation between modulating voltage and modulation in- 
dex, ft i ' ^ en f rom Eq. (7-16), 




(7-16) 



smaller value of Eg, which in turn tends to make £3, increase. 
Specifically, in the present system, ^ £ / , and for S^= 10, 





(7-17) 
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Also, from Eq. (7-16), for the conditions given above 

Pr = Pi? Pi ~ — m 

Thus, k^ = 4.9, where C is the value of E<> corresponding to 
C 

S-j_ = 10, keeping in mind that P^.= constants 
Now from Eq. (7-15), solving for Eg 






S ! ~t~ / 



(7-18) 



and, 



C={^X^)-'\ It? , 

from which 

• 

Substituting this quantity for k-^ and using Eq. (7-18), Eq. 
(7-17) becomes 



P,= 



i+- v 



’{£-i 



S,+/ 



/ + dT/.V 



i/il 



(7-19) 



Evaluating Eq. (7-19) for S^= l/lOO or -20 decibels, 

& = L = = ./63 . 

7 + SL 

Prom Eq. (6-1) it is noted that for large R, increasing ^3 
from .02 to .163 improves the signal-to-noise ratio S 2 by 
approximately — ^><6- 5" or IQ.3db,. 

It remains to determine the effect of the I.P. amplifier 
phase shift characteristic upon the phase of the modulating 
wave at the output of the discriminator. Disregarding the 
noise voltage, the wave at the input to the first I.P. stage 
may be written (12) 
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^4 /4 -h A^n- oJ ( 

CO 

~ A ■+ nc^)^] 



*--/ 

oo 



+ ft J" 'i ( pi) - «<^, ) t| | 

l~t 



( 7 - 20 ) 



r\--l 



The phase shift and amplitude response vs. frequency may 
be obtained from Eq. (3-8) p. 42 of Terman (22), 



© = , 
('/•y-rf ' 



\ 



(7-21) 



Mfl _ 


/ 


K| 


/y 6 y (/+S) 



where 



S=T- 

* n 



Then the wave at the output of the I.F. amplifier may be 
written 

e o (i) =£,jZ(^^n.(^t-e') 

OO 

-h V6 0 J n (p^^£cJ 0 +„ +■ 4> t J 



n-t 

co 



+ 



y b n (-0 4-<f> n jl (7-22) 

r£/ ^ 
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where b n , <^> and b^, <p^ are the amplitude responses and 
phase shifts at the frequencies f Q + nf^ and f Q - nf-^, respec- 
tively. 

Values of are tabulated below for = .02 and 

.163 and n = 0, 1,2,3* All higher orders are zero through at 
least the fifth decimal. 



Jn\f' 


CM 

O 

• 

l 


.163 


J 0 




799990 


' .99337 




.01000 


.08123 


’ J 2 


~ .00605 


.00332 


[ 

\ 

[ 


.00000 


.00001 



Values of b n , (f> n , b^ and may be obtained from Eq. (7-21) 
and are tabulated below. 



n 


0 


1 




b 

n 


1 

I \ 


.843 




.609 | *400 


K 1 


.970 

L . . - i 


.852 | .606 

. . . .. i 


r 

! 

C* 

"Ok 

i : 


0 ! 57.0° 1 109* 0° ; 154.1° 

1 1 . : . 


i & ! 0 


-60.0° : -118.9° • -177.3° 

1 } 



Taking p = .02, it is seen from Eq. (7-22) that only the 
first order sidebands are significant in representing frequency 
modulation. A phasor diagram of this situation is shown below 
in Fig. (7-3)* 
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Fig. (7-3) 

It can be seen from the diagram and Eq. (7-22) that the 
phase of the I.F. wave will pass through its zero position 
when the two first order sidebands are in such a position that 
their phasor sum has no component perpendicular to the u5 £ Hb 
phasor. This condition is expressed mathematically by writ- 
ing 

.8v3^(?.^V*)= (7-23) 

from which 

&orL f/' $ -/■*) - /. I S’ Atsn* 

C -&& X 5 Ak X — (eiS~ Ak^/,5 C&A. X — /< /-S £" At'?'' X 

. /i " /, x - <. /b" /. £ Au,-,^ X 

/ • & 

J2Ltx.X - — — mU-v IS 

7.,/t 

- ,ooi?i 

X k (oS . (7-24) 
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This is the amount by which the phase of the wave is shifted 
from the position it would have if the amplitude response were 
perfectly symmetrical about f Q . Even if it were not negli- 
gible, it could be calibrated out. It is recorded here simply 
to illustrate the method of calculation and to serve as a com- 
parison with the phase shifts under different conditions. 

Again it is assumed that the I.F. amplifier is correctly 
tuned to f 0 and is increased to .163. The sidebands at 
f Q — 2f^ must now be considered. A phasor diagram of this 
situation is shown below in Fig. (7-4). 




Fig. (7-4) 

It appears intuitively that the effective phase shift con- 
tributed by the unsymmetrical second order sidebands will be 
negligible because of their small amplitudes. To show that 
this is true, assume that the smaller one is removed 
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completely. The first order sidebands would have a nullify- 



value of x less than 10 of arc different from that of Eq. 
(7-24). Since the system accuracy as determined by the servo 



phasemeter is no better than 1.5 feet and 10 of arc corres- 
ponds to approximately 0.92 feet (at the ranging frequency of 
491 Kc) , the shift in zero axis crossing due to the second 
order sidebands is obviously negligible. 

The above investigation assumed the I.F. amplifier to be 
correctly tuned to a center frequency f Q . Now it will be 
assumed that the carrier frequency remains at f Q but that the 
center frequency of the I.F. amplifier is shifted upward to 
f 0 +f^. This is an arbitrary choice but is based upon the 
fact that the 3db bandwidth is approximately 4fi> hence a 
shift in amplifier center frequency of f^ marks about the 
maximum shift allowed for in the design "safety factor". The 
amplifier was of course designed for {3 t - .02, so that only 
the first pair of sidebands had to be considered. 

Taking @ = .02 and again disregarding the sidebands at 
f 0 ± 2f^, Eq. (7-22) may be written 



ing component perpendicular to the reference for a 



Q 0 (£) - 3j.9j J 
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Using the method of Eqs. (7-23) and (7-24) it is found that 



X - 




(7-26) 



Now when £ - .163 it is seen that the amplitudes of the 
second order sidebands stand in almost exactly the same ratio 
as in the first case considered, while their phase angles 
differ by only 0.3° • Thus, the effect upon zero crossings of 
the second order sidebands is again negligible. 

The amplitude modulation which is present when |9 = .163 
is at twice the modulating frequency f^ and is not passed by 
the tuned circuits of the compensation amplifier. 
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8. Conclusions. 

While the analysis of Section 5 is far from rigorous, it 
is considered a fair conclusion that use of the phase discrimi- 
nator in place of the frequency discriminator would provide a 
significant extension of system range. The experimental re- 
sults of Section 6 support this conclusion. In fact, the 
experimental results are likely to prove more typical in other 
discriminators of similar construction, since the assumptions 
and approximations involved in the work of Section 5 are 
really quite severe. The assumption that the diodes are per- 
fect switches, for example, is not realistic when the signal 
is very small. Linear operation of the diodes is similarly 
unlikely when the signal is very large. Complete symmetry of 
components is still another condition which is not attainable. 
However, it is believed that the mathematical model is a 
reasonable approximation in the region where the r.m.s. refer- 
ence voltage is maintained very large compared with the total 
r.m.s. voltage in the I.F. amplifier [conditions of Eq. (6-2)] 
or carrier-to-noise ratio (S-^) is large. This condition is 
quite easily met when total r.m.s. voltage in the I.F. ampli- 
flier is maintained nearly constant. The reference voltage 
amplitude is simply kept constant at a value which is large 
with respect to the carrier amplitude where S-^ = +10 db. This 
is a strong point in favor of using the ’’constant total power” 
AGC method of Section 7. 

The desirability of keeping circuit changes to a minimum 
suggests another advantage of retaining the present AGC 



45 



system rather than using a "carrier-only" type. Even the ad- 
dition of an amplifier for the AGO voltage would be relatively 
simple compared with the extra circuitry required, for ex- 
ample, in Pig. (1-1). 
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9. Recommendations for Further Work. 

It appears that the phase discriminator will find many 
applications in this and other areas of research where phase 
measurements and/or correlation detection are employed. Radio 
navigation, radar, color television and many other electronic 
devices could be cited as examples. Consequently, a really 
thorough investigation of the phase discriminator might be a 
highly useful project. 

It was mentioned in Section 4 that many non-linear de- 
vices have been investigated theoretically to determine their 
effects on the spectrum of signal plus noise. Regardless of 
the method used, a fairly common occurrence is the appearance 
of an integral or infinite series which cannot be expressed 
in closed form. 

Blachman (2), Middleton (17) » Rice (19) and Lawson and 
Uhlenbeck (16) are a few who have reported such results. 

Their procedure is then to assume "very large carrier-to-noise 
ratios" and "very small carrier-to-noise ratios" which reduce 
the equations to closed forms. These results are useful for 
some practical applications such as the detection of extremely 
weak signals in noise or circuit response to signal alone. 
However, equally valid results are often obtainable by the 
methods of modulation theory for these limiting cases. The 
region in which calculations are most difficult (carrier-to- 
noise ratios roughly -20 db to +20 db) would therefore be the 
region of most practical interest. 

It would then seem that a logical approach to a thorough 
investigation of the phase discriminator problem would involve 
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a determination of an integral or series form most amenable 
to solution on a digital computer. Two series might be re- 
quired to cover the entire range (17) but as long as they 
could be derived from the initial model without approximations 
the numerical results could be obtained to any degree of ac- 
curacy desired. 

Eq. (4-1) requires the knowledge of F( ju) for the device 
being considered. Rice (19) has tabulated F(ju) for a number 
of non-linear devices in Appendix 4A of his paper, but no one 
has written this function for the phase discriminator. This 
function could probably be calculated with the aid of the 
papers by Rice (19) » Rice and Bennett (20) and Middleton (17). 
Having obtained this F(ju), the simplest form for numerical 
analysis would have to be reasoned out. 
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APPENDIX A 



A . 1 AVERAGE VALUES 



The characteristic function for a random variable y, or 
for the corresponding distribution, is defined by (10) 

CK> 






d y- W t (y) £ 



j 



"t 






(A. 1-1) 



'-cc 



This result may be generalized for a distribution of n 
variables (10) but for the present purpose it is sufficient 
to write 




Since t^ and t^ are constants as far as the integration 
is concerned, it is possible to differentiate both sides with 
respect to t^ and t^. Differentiating m times with t ( and n 
times with t the result is 




If now t^ and t 

sT*" 'pQt, , A) 

d"Z 2 



are set eaual to zero, the result is 
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The right side of this equation is just (^J 
Hence , 



r r = 



-hi i-o 

/J\ 


itY) 4- / / / / N -i 

£ 


W 


— / ^ J 



ir t -o 

£ k --o , 



(A. 1-3) 



A. 2 CHARACTERISTIC FUNCTION FOR W, 



The characteristic function for will now be calcula- 



ted . 



oo 



'/'Ov j l Z 7T - 















4 J h 
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ZG^O-f) 



Jdj t 



The integration over is performed first. Completing the 
square in the exponent the result is 

OO 






/-CO 






oc< 



2cr^(/-^y 



/- Co 



50 



Now making the substitution 



Z-rrz = 2 , ' 

this integral reduces to 

2 7T- **r [-£$=?} +S*, 



p 2 






OO 




dz 



c -*r[i 



cTJ 



■ H7T Z Z Z 

^(hf) 



OO 



Making use of pair No. 708.0 Campbell and Foster (6) this be- 
comes 

- Z Z . 2\ 

a 



^7T <r z (/-f z ) - exp 



,2. ,2. 



4- j 



and 



oo 






fw^ Jd**- 6xp d^~d*o~ f) 




2 - 

2 



M&z+pi.) - jr-i_ 



€ 



Setting 



= 2 . 7 r£ 



and again using pair No. 708.0, this becomes 
Vz7 rcr 2 - 

which simplifies to 

- ^\±f'+Zoi l t- z +±^ 

f . (A. 2-1) 
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APPENDIX B 



Transformation of Eq. (5-7) requires the determination 
of the Fourier transforms of C<£<E uj ( 'V ^ p C&Xs cO f 'Y 

, f ( Y end o 3 - 



As a starting point, it is known that (16) 
00 

d± 6 




-jzjr-ft 




(B-l) 



-OO 



where 6(f) is the Dirac delta function and has the following 
properties : 



CO 




^00 



- / ( 

<5W = ctf-f) 



So 



7^= 2 






To transform Gxi a), X , it is first written 



C6^u\X = — 
£ 



So 



3* 



r^ = f 

C/..C 






-jdJ. v / 

cj, r 6 T 



00 






CO 
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_L 

2 



6(+-f, 



-f 



2 



<5(4 + -Q 



or 



S 7 * S(-f~-f,). (positive frequencies) (B-2) 
The transform of p(T) is obviously 



cr 



(B-3) 



O 



The transform of p (T) is most easily taken by using 
the convolution theorem (25)* Thus, 

0 ° A . / oo 

T Jr " -pt jGJlUf-jdA 



=. A-n-£l 

bL bJ 



(B-4) 



o<-f < 8 



Denoting the transform of p 2 (T) by and using 
the convolution theorem again, the transform of p 3 (T) is 
found to be 




CT' 



oo 

f H(^) (f-2) cIa 



/-oo 



O < + < B /x_ 
ss/Sci m h <3. r , 



(B-5) 
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It was seen above (l.q. B-2) that transforming CsrO- (*o K X 
shifted the spectrum of the unit impulse by - -f . In an 
entirely similar manner it is found that 




-jvrtr 






dr 










) 




(B-6) 




o < -f <.-f t 

■f, < -f < 6-1 ; 

6-f,< (<&+{, 



( B-7) 
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APPENDIX C 



Values of I (x) and I n (x) for O 4 X+. 5" and of £ -XA) 
o 1 

and € for ^ x ^-ZO may be obtained from British 

Association for the Advancement of Science Mathematical 
Tables, Volume VI . (5). 

The expression 




£ * 

V 7 ^ 



(i -f-^Iod) + E I, ( 




was evaluated in 3 db increments of S^ from -15 db to +15 
db. Results are tabulated below: 



S 1 (db) 


s i 


T 


-15 


.03125 


1.000 


-12 


.0625 


1.000 


-9 


.125 


1.000 


-6 


.25 


1.000 


-3 


.5 


1.008 


0 


1 


1.024 


3 


2 


1.045 


6 


4 


1.073 


9 


8 


1.099 


12 


16 


1.110 


15 


32 


1.120 


The asymptotic 


values of I Q (x) 


and I (x) are found on 
1 


p. 271 of the same 


tables. Substituting these expressions 


into the right side 


of the equation 


for T(S ) and letting 
1 
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= X , the result is 
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Now letting X 



a?and evaluating the limit, 
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